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Abstract 

We demonstrate saturable absorber behavior of n-type semiconductors GaAs, GaP and Ge in 
THz frequency range at room temperature using nonlinear THz spectroscopy. The saturation 
mechanism is based on a decrease in electron conductivity of semiconductors at high electron 
momentum states, due to conduction band nonparabolicity and scattering into satellite valleys in 
strong THz fields. Saturable absorber parameters, such as linear and non-saturable transmission, 
and saturation fluence, are extracted by fits to a classic saturable absorber model. Further, we 
observe THz pulse shortening, and an increase of the group refractive index of the samples at 
higher THz pulse peak fields. 
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Semiconductor saturable absorbers and saturable absorber mirrors (SESAMs) are rou- 
tinely used for ultrafast laser modelocking and ultrafast signal control [1]. Saturable ab- 
sorbers operating in the visible and infrared wavelength ranges rely on bleaching of two-level 
electronic systems, usually realized by an interband transition in semiconductor quantum 
wells or quantum dots i.2i]. Clearly, such quantum- confined semiconductor systems will not 
be suitable for applications in the far-infrared (THz) range, where the photon energy is much 
smaller than the bandgap energy of most semiconductors, and where thermal population of 
such closely spaced electronic levels would dominate. The main loss factor in doped semicon- 
ductors in the THz frequency range is free-carrier absorption. The THz signal attenuation 
is roughly proportional to the conductivity of the material a = | e | [/ie'^-e + /^/i^/i], where e is 
elementary charge, and /ie,/i and n^^h are the mobilities and concentrations of electrons and 
holes, respectively js, 4|. 

In this Letter we demonstrate saturable absorbers for the THz frequency range, based on 
n-type bulk semiconductors, where the carrier mobility is modulated by nonlinear electron 
transport caused by the THz electric field, thus affecting the conductivity of the sample. 
Normally, the application of external electric fields leads to acceleration of carriers in the 
lowest-energy valley of the conduction band of an n-type semiconductor. At high-momentum 
states the valley nonparabolicity becomes pronounced, which leads to an increase in the 
effective mass and thus to a reduction of the mobility Ue, consequently leading to a reduced 
d*c..o 10. in *e TH. .a„,e g M e„ou,K e.ect.c fields, l^ervalley is 
possible 0, 0] , leading to electron transfer into the satellite valleys with reduced curvature. 
This again results in a smaller electron mobility compared to that of the conduction band 



minimum, and consequently lower THz dielectric loss. Detailed studies of time-reso 



high-THz-field transport in bulk semiconductors were recently published in Refs. |8NlO|. 
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In our experiment, we generated high-power single-cycle THz pulses by tilted pulse- front 
optical rectification in a lithium niobate crystal of 800-nm, 80-fs laser pulses provided by a 
1-kHz repetition rate Ti:Sapphire amplifier Uj. The THz pulses were coUimated and then 
refocused onto a sample point using a set of off-axis paraboloidal mirrors. After propagation 
through the sample point, the THz pulses were guided to a 0.5-mm thick undoped <110>- 
oriented ZnTe crystal for detection using standard free-space electro-optic sampling (FEOS) 
[l^ . A pair of wire-grid polarizers was introduced into the THz beam path before the 
sample point, which allowed us to controllably attenuate the THz signal by adjusting the 
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FIG. 1: (a) Time-resolved electric field of a THz pulse with a fluence of 50 ^J/cm^, and a peak 
electric field of 300 kV/cm. Corresponding amplitude spectrum with noise floor indicated, (b) 
Instantaneous THz intensity at the sample position calculated from the square of the measured 
electric field. 

angle between the polarizer axes. The maximum THz pulse energy was measured to be 1.5 
fiJ using a calibrated pyroelectric detector, in the THz beam waist of 1.0 mm (intensity, 
1/e^) at the sample location. This is equivalent to a maximal THz fluence of Fmax = 50 
/iJ/cm^, corresponding to a peak electric field of 300 kV/cm. The peak field was determined 
by calibrating the squared time-domain electro-optic signal with the known value of pulse 
fluence. Fig. [1] shows the temporal dependency of electric field as measured by FEOS, the 
amplitude frequency spectrum, and corresponding instantaneous intensity at the position 
of the sample of the strongest THz signal used in our experiments. The samples under 
study were n-type semiconductors: GaAs with a carrier concentration of 8 x 10^^ cm~^ 
and a thickness of 0.4 mm; GaP with a carrier concentration of 10^^ cm~'^ and thickness 
of 0.3 mm; and Ge with a carrier concentration of 10^^ cm~^ and a thickness of 6 mm. 
The measurements were carried out at room temperature. Fig. [2] shows the field and 
power transmission coefficients as a function of THz pump field and fluence respectively 
in GaAs, GaP and Ge. The field and power transmission coefficients were obtained by 



3 



0.32 
0.28 
0.24 
0.20 
0.16 



c 




a 








E 

0) 


0.084 


o 




u 




c 


0.080 


o 




'ui 




in 




1 


0.076 


in 




c 




IS 

1_ 


0.072 


♦J 




■a 






0.48 


E 






0.46 




0.44 




0.42 



a) 


• 




GaAs . • 

• 




• 

■ 




, • * d = 0.4 mm 




rig = 8 10 cm"" 


c) 


• 




GaP 

• 

■ 




• 

^ ^ • d = 0.3 mm 


• 


rig = 10^^ cm-^ 


e) 






Ge . . • • * 

• 

• 


% 


• 

• d = 6 mm 


•. 


ng = 10^''cm-3 



50 100150 200 250 300 

Peak THz field [kV/cm] 




0.10 
0.08 
0.06 



T„ = 0.148 
T|,„= 0.024 

F„,= 8.2 nJ/cm2 0.02 




T„, = 0.23 
T|i„ = 0.16 

F.„= 3.1 nJ/cm2j0.17 
10 20 30 40 50 

THz pulse fluence [nJ/cm^] 



FIG. 2: Symbols: field transmission coefficient as a function of peak THz pulse field, and power 
transmission coefficient as a function of THz pulse fluence for GaAs (a-b), GaP (c-d), and Ge (e-f). 
Solid lines - saturable transmission function fit to the measured power transmission coefficients. 
See text for details. 



integrating either the modulus or the square, respectively, of the transmitted THz fields 
transmitted through the sample, and dividing them by reference values recorded without 
the sample in the beam path. In all our samples we observed increased transmission at higher 
pump fiuences. In particular, we observed a nearly five-fold increase in power transmission 
coefficient for GaAs sample in the THz pulse fiuence range used in our experiments. The 
sohd lines in Figs. |3]|^b,d,f) are fits using a saturable power transmission function, defined 
after Ref. [l^ as 

TiFp) = T^s— ^^V^^ (1) 

where Tu^ and T^s are linear and non-saturable power transmission coefficients, Fp is the 
pump fiuence, and Fsat is the saturation fiuence. Using fits with Eq. [T] we were able to 
extract the saturable absorber parameters for our semiconductor samples. These parameters, 
namely linear and non-saturable transmission, and saturation fiuence are indicated in Fig. 
|2l In particular, the saturation fiuence F^at was found to be 8.2 /iJ/cm^ for GaAs, 20.9 
yuJ/cm^ for GaP, and 3.1 yuJ/cm^ for Ge. Interestingly, at the lowest pump fiuences the 
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transmission function for Ge actually decreases slightly, thus demonstrating the opposite, 
optical limiting behavior. This can be possibly explained by the fact that in Ge the high- 
mobility F-valley is not initially populated, and thus initial sample conductivity is lower, 
as discussed in Ref. fsl- In all our samples we observed shortening of the THz pulses, as 
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FIG. 3: (a) Shape of the THz pulses before (reference) and after (sample) propagating through the 
0.4-mm thick GaAs sample. This reference pulse had peak electric field strength of 292 kV/cm. 
(b) Modulus of Hilbert transforms of these THz pulses. Horizontal arrows indicate FWHMs of 
the pulses, and vertical arrows indicate the mean-weighted arrival times of the pulses, (c) Pulse 
shortening factors and (d) group refractive indices of GaAs, GaP, and Ge samples as functions of 
peak electric field of incident THz pulse. See text for details. 

they propagated through the samples, with increase in peak electric field of the THz pulse. 
Such pulse-shortening behavior is also characteristic for saturable absorbers in the optical 
range. In order to quantify the shortening of single-cycle THz pulses used in our experiment, 
we calculated a pulse shortening factor defined as the ratio of full width at half maximum 
(FWHM) of modulus of Hilbert transforms of sample and reference pulses, as illustrated 
in Fig. [3]^a,b). The dependency of the pulse shortening factor on the peak electric field 
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of the incident THz pulse is shown in Fig. [3]^c). For GaP and Ge the maximum relative 
pulse shortening (i.e. the difference between largest and smallest shortening factors) reaches 
approximately 5%, and its dependency on THz peak field is more or less linear. For GaAs, 
however, this dependency appears to have a threshold at THz peak field of around 150 
kV/cm and, and the relative pulse shortening reaches the value of approximately 10%. 

Saturable absorption and pulse shortening are accompanied by an increase in group re- 
fractive index in all three samples, as shown in Fig. |3t^d). The group index Ug at various 
pump intensities was calculated using the difference of arrival times Ar between sample 
and reference pulses and taking into account the sample thickness d by using the relation 
Hg = ATc/d+ 1. The arrival times were obtained from mean-weighted maxima of the mod- 
ulus of THz pulse Hilbert transforms. We observe a clear saturating growth trend in the 
group refractive indices of all our samples with increasing THz fields. The most dramatic 
index growth is observed in GaAs, and it reaches Aug = 0.14, whereas for Ge it is only 
Ang = 0.002 in the full THz pulse peak field range available in our experiments. We note 
that very similar dependencies were obtained if the maxima of THz pulses, and not of the 
modula of their Hilbert transforms were analyzed. 

The observed increase in group index with stronger THz pump is indeed anticipated, 
and two contributions to this effect are considered. The first contribution is the change in 
reflectivity at the sample surface {s]. A decrease in the sample absorption, resulting from 
the strong THz excitation, leads to a positive phase shift of THz pulse at the interface 
between the air and absorbing sample. This causes a temporal delay of the transmitted 
THz pulse and consequently a larger group refractive index. The other contribution to a 
group index growth with stronger pump stems from the propagation through the bulk of the 
sample with the phase index dispersion in the frequency range below resonance (in this case 
- plasma resonance). Larger effective mass m*, acquired by the carriers at high-momentum 
states as a result of intense THz excitation, will lead to a decrease in plasma frequency 
i^p = (jT-eC^/eom*)^/^, thus making the frequency dependency of phase index n^{uj) more 
steep. This would result in the higher value of a group index = n,^ + udn,j)/duj for 
the stronger excited samples. The same effect will arise from the effective reduction of the 
density of highly mobile carriers rif. at the bottom of the F- valley. Both interface and bulk 
contributions to group index growth obviously depend on the strength of the THz excitation 
and will thus saturate with the overall THz transmission, which is in accordance with our 
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observations. We note that a similar increase in the pulse delay and thus in the group index 
with stronger optical pump was also observed in saturable absorbers for the optical range, 



as described in Ref. 



1^. 



We further note that n-type GaP is often used as a broadband FEOS detector, where 
parasitic etalon reflections of the THz pulses within the crystal are conveniently suppressed 
by the free-carrier absorption. As we have shown here, in the case of high THz field strengths 
the results of THz spectroscopy will be distorted due to n-GaP nonlinearities. Therefore, for 
etalon suppression in high-field experiments it is advisable to use combined active-inactive 
crystals [15] instead of doped crystals. 

In conclusion, we have observed ultrafast saturable absorption in n-type semiconductors 
GaAs, GaP, and Ge in the THz range at room temperature. We were able to extract linear 
and non-saturable transmission coefficients, as well as the saturation fiuence using a standard 
saturable absorber model developed for the optical frequency range. Our data shows THz 
pulse shortening, and growth of group refractive index of our samples with increase in 
THz excitation, correlated with the saturating transmission through the samples. Although 
the origin of the absorption mechanism in the THz and in the optical frequency ranges is 
completely different, we note that the saturation fluences observed here are within the same 
order of magnitude (i.e. few to tens of /iJ/cm^), as the values reported for SESAMs in the 



optical range, such as e.g. quantum dot SESAMs [2, ll6|. Optimized THz SESAMs could 
potentially be used with the THz sources like quantum cascade lasers (QCLs), leading for 
example to external-cavity mode-locked QCLs delivering high-energy ultrashort THz pulses. 
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Figure captions 

Caption to Fig. 1 

(a) Time- resolved electric field of a THz pulse with a fiuence of 50 /iJ/cm^, and a peak 
electric field of 300 kV/cm. Corresponding amplitude spectrum with noise floor indicated. 

(b) Instantaneous THz intensity at the sample position calculated from the square of the 
measured electric field. 

Caption to Fig. 2 

Symbols: field transmission coefficient as a function of peak THz pulse field, and power 

transmission coefficient as a function of THz pulse flucncc for GaAs (a-b), GaP (c-d), 
and Ge (e-f). Solid lines - saturable transmission function fit to the measured power 
transmission coefficients. See text for details. 
Caption to Fig. 3 

(a) Shape of the THz pulses before (reference) and after (sample) propagating through the 
0.4-mm thick GaAs sample. This reference pulse had peak electric field strength of 292 
kV/cm. (b) Modulus of Hilbert transforms of these THz pulses. Horizontal arrows indicate 
FWHMs of the pulses, and vertical arrows indicate the mean-weighted arrival times of the 
pulses, (c) Pulse shortening factors and (d) group refractive indices of GaAs, GaP, and Ge 
samples as functions of peak electric field of incident THz pulse. See text for details. 
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